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Abstract 

We present a N = 1 Supersymmetric extension of a spin-charge separated effective SU(2) x Us(l) 
'particle-hole' gauge theory of excitations about the nodes of the gap of a d-wave planar magnetic 
superconductor. The supersymmetry is achieved without introducing extra degrees of freedom, as 
compared to the non-supersymmetric models. The only exception, the introduction of gaugino fieds, 
finds a natural physical interpretation as describing interlayer coupling in the statistical model. The 
low-energy continuum theory is described by a relativistic (2+l)-dimensional supersymmetric CP 1 
(T-model with Gross-Neveu-Thirring-type four-fermion interactions. We emphasize the crucial role of 
the CP 1 constraint in inducing a non-trivial dynamical mass generation for fermions (and thus super- 
conductivity), in a way compatible with manifest N = 1 supersymmetry. We also give a preliminary 
discussion of non-perturbative effects. We argue that supersymmetry suppresses the dangerous for 
superconductivity instanton contributions to the mass of the perturbatively massless gauge boson of 
the unbroken U(l) subgroup of SU(2). Finally, we point out the possibility of applying these ideas 
to effective gauge models of spin-charge separation in one-space dimensional superconducting chains 
of holons, which, for example, have recently been claimed to be important in the stripe phase of 
underdoped cuprates. 
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1 Introduction 



In ref. it was argued that the doped large-?/ Hubbard (antiferromagnetic) models possess 
a hidden local non-Abelian 577(2) x Us{l) phase symmetry related to spin interactions. This 
symmetry was discovered using an appropriate 'particle-hole symmetric formalism' for the elec- 
tron operators |2[], and employing a generalized slave-] ermion ansatz for spin-charge separation, 
which allows intersublattice hopping for holons, viewed as fermions. The spin-charge separa- 
tion may be physically interpreted as implying an effective 'substructure' of the electrons due 
to the many body interactions in the medium. This sort of idea, originating from Anderson's 
RVB theory of spinons and holons H, was also pursued recently by Laughlin, although from a 
(formally at least) different perspective J3J. 

In ref. [TJ we have argued in favour of the opening of a fermion gap at the nodes of 
a d-wave gap of a superconducting antiferromagnet. Linearization of the fermion spectrum 
about such nodes leads to a relativistic Dirac spectrum for holons, with the role of the limiting 
velocity being played by the fermi velocity 0, [L|. Such systems might be of relevance to the 
physics of high-temperature superconductors, since recently it is believed that high-temperature 
superconductivity in cuprates is highly anisotropic and the gap symmetry is of <i-wave type ||, 
with the gap vanishing along lines of nodes on the Fermi surface []. 

The key suggestion in ref. |TJ, which lead to the non-abelian gauge symmetry structure for 
the doped antiferromagnet, with the constraint of not more than one electron per lattice site, 
was the slave-fermion spin-charge separation ansatz for physical electron operators JTJ: 

_ ( Ci C 2 \ _ ? ^ ( Ipl lp 2 \ ( Z X -Z 2 \ m 

*<* = U -c\)=^'^={-4 a) t { * * J, < l > 

where % is a lattice site index, c a , a = 1,2 are electron annihilation operators, the Grassmann 
variables ipi, i = 1,2 play the role of holon excitations, while the bosonic fields Zi,i = 1,2, 
represent magnon (bosonized spinon) excitations ||. The ansatz (p]) has spin-electric-charge 
separation, since only the fields ipi carry electric charge. 

As argued in ref. |IJ the ansatz is characterised by the following hidden local phase (gauge) 
symmetry structure: 

G = SU{2) x U S {1) x U E {1) (2) 

The gauge SU(2) symmetry pertains to the spin degrees of freedom. The local Us(l) 'statistical' 
phase symmetry allows fractional statistics of the spin and charge excitations. This is an 
exclusive feature of the three dimensional geometry, and is similar in spirit to the bosonization 
technique of the spin-charge separation ansatz of ref. ||11|| . Finally the ?7e(1) symmetry is due 
to the electric charge of the holons. 

1 There is also recent experimental evidence on the possibility of the opening of a gap at such nodes, triggered 
by either magnetic fields [Q or by magnetic impurities Q , and although such phenomena might admit alternative 
(more conventional) explanations |l^| , however the role of spin-charge separation in this context still remains 
a challenging project. 
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It is the purpose of this work to discuss the possibility of a hidden supersymmetry in the 
ansatz ([!]). Note that supersymmetric extensions of J = ±2t models for doped antiferromag- 



nets, in one and two spatial dimensions, have already appeared in the existing literature |L2 
even in the context of spin-charge separated anyon models [T3|. However, as far as we are aware, 
such supersymmetries have not been associated so far with any specific dynamical properties 
of the antiferromagnet. In contrast, in our approach here, based on the non-trivial ansatz (|I|), 
the supersymmetry constitutes a non-trivial dynamical property of the spin-charge separated 
vacuum for holons and spinons, by viewing them as supersymmetric partners. Due to the rich 
group structure (0), many possibilities arise in the study of the phase diagrams of these theories, 



in the context of the modern perspective advocated in the work of Seiberg and Witten [14, 15 



In particular, duality symmetries in the infrared region of the supersymmetric model, connect- 
ing various theories with the same non-trivial infrared fixed-point |15| , may prove very useful in 
a renormalization-group study of the dynamics of the gauge fields in both, the superconducting 



and the normal phases of the model, in the spirit of ref. [Jig] . The important issue is that the 
introduction of N = 1 supesymmetry, hidden in the spin-charge separation ansatz (|l|), does 
not require the introduction of unphysical degrees of freedom. As we shall see, the only extra 
degrees of freedom, as compared to the non-supersymmetric case [p]] , are the gauginos of the 
local hidden gauge symmetry, which, however, admit the natural interpretation of describing 
interlayer hopping of spin and charge degrees of freedom (hopping of 'real' electrons). 

We should stress that, within a condensed- matter context, the supersymmetry refers to 
the relativistic field theory at the nodes of a ci-wave superconducting gap 0. In this sense, the 
supersymmetric dynamics of the spinons and holons would require equality of the spin gap with 
the fermion (superconducting) gap at such nodes. At a microscopic level, this would imply some 
particular relation among the microscopic parameters of the model, such as hopping matrix 
elements and Heisenberg interactions. This calls for comparison with the J = ±2t special point, 
where the graded (supersymmetric) algebra in the spectrum of the doped antiferromagnets 
appears [ |T2"| , [HJ. However, as we shall see, the situation in our case is more complicated, since 
there are more parameters entering the dynamical scenario of the gauge theory based on the 
spin-charge separation ansatz (pi). 



2 Review of the (continuum) model and its supercon- 
ducting properties 

Before embarking to a description of the supersymmetric extension we consider it as useful to 
review first the properties of the statistical model of ref. [0], some of which will be crucial 
for the supersymmetric extension. The pertinent long-wavelength gauge model, describing the 
low-energy dynamics of the large-U Hubbard antiferromagnet in the spin-charge separation 
phase ([I]), can be cast in a conventional relativistic lattice gauge-theory, provided one changes 

2 Galilean supersymmetry, as symmetry of the spectrum between bosonic and fermionic degrees of freedom, 
may also occur away from the nodes. This is left for future work. 
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representation of the SU(2) group, and, instead of working with 2x2 matrices, one uses a 
representation in which the fermionic matrices ip a p are represented as two-component (Dirac) 
spinors in 'colour' space: 

= (Vi - Va) f > *2,< = (*h . , i = Lattice site (3) 

By assuming a background field of flux tt per lattice plaquette ||, and considering 

quantum fluctuations around this background for the Us(l) gauge field, one can obtain the 
conventional lattice Dirac action for the fermion excitations about a node in the fermi surface J5], 

01 

In the above context, a strongly coupled Us(l) group can dynamically generate a mass 
gap in the holon spectrum |T^, [5], |T^, |2D|, |2T|, which breaks the 577(2) local symmetry down 
to its Abelian subgroup generated by the a 3 Pauli matrix [I], p2[ . From the view point of 
the statistical model of ref. Q, the breaking of the SU(2) symmetry may be interpreted as 
restricting the holon hopping effectively to a single sublattice, since the intersublattice hopping 
is suppressed by the mass of the gauge bosons. 

The (naive) continuum limit of the low-energy theory about such nodes on the fermi surface 
of the planar antiferromagnet, then, is described by a CP 1 model coupled to Dirac fermions || 

C 2 = gt\(d, - (g 2 /gi)BX - a,)z\ 2 + WDtf^ (4) 

where now = <9 M — ia^ — %{gij g\)o~ a B aiL — - C A^ gf, i = 1, 2 have dimensions of mass, B° is 
the gauge potential of the local ('spin') SU(2) group, generated (in two-component notation 
for fermions ) by the Pauli matrices a a , the £7s(l) ('fractional statistics') field, and is 
an external electromagnetic potential, which will be ignored in the subsequent discussion. In 
terms of the microscopic model, g\ ~ J5, where J is the Heisenberg exchange energy, and S 
is the doping concentration. An important ingredient in the above formalism is the no-double 
occupancy constraint, which in terms of the z and ty a , a = 1,2, fields, with a a 'colour' index, 
can be written as: 

y £[3 a z a + 0* a a 3 * a ] = l (5) 

a=l 

where 03 acts in spinor space, and the fermions \& are viewed as two-component spinors, related 
to the spinors \I/ (|3p by appropriate rescalings so as to ensure the canonical kinetic (Dirac) 
term 0. This results in the presence of the constant (3 (with dimensions of [mass] -2 ) in the 
constraint (|5|) |nj. In the context of the microscopic model, these constants are expressed in 
terms of the hopping and Heisenberg exchange energies [0, |TJ, and one has that \/3\ << 1. It 



can be shown |[] that the constraint (|5|) is essential in ensuring the consistency of the ansatz 
([!]) with the canonical commutation relations of the electron operators. 

The presence of the ty'ty (non-relativistic) fermion number term in the constraint (|5|) appears 
at first sight to complicate things, since the conventional CP 1 constraint \z\ 2 = 1 is no longer 

3 In the model of ref. Q, due to the Dirac nature of the resulting spinors, and \& are viewed as independent 
variables in a path integral, which implies that one can redefine ty> — > 
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valid. However, these extra terms can be rendered inocuous for the dynamics of the effective 
theory. Indeed, by integrating out the (non-propagating ) gauge fields in (f|) we obtain |23| : 



|-[1 - (3^ a a^ a ]- l Ti {z a d^z a - z a d^z a - ig^l^l + 0# a ) 2 + 

61n[l - pW(T 3 * a ] (6) 

where the last term is absent in the usual CP 1 models. Expanding this term in powers of the 
(small) parameter /3 « 1, one obtains: 

61n[l - 0Ta 3 ^ a ) ~ -60 a a 3 V a + 3/3 2 (FV 3 * a ) 2 + . . . (7) 

where the . . . indicate six- and higher order -fermion contact terms, not renormalizable, even 
in large-N limits, which constitute irrelevant operators, in a renormalizat ion-group sense, not 
affecting the low-energy (infrared) structure of the effective theory, we are interested in. 

Applying a Hartree-Fock linearization to the four-fermion interactions, one obtains terms 
of the form: 

3(3 2 < ^03*0 > ¥V 3 * Q (8) 

Collecting the Wty terms together, one then obtains a fermion- density term in the effective 
lagrangian of the form: 

Lfx = (-6/3 + 3/3 2 < FV 3 * Q >)^V 3 ^ (9) 

Upon inserting the constraint (g) via a Lagrange multiplier field X(x) in the path integral, one 
may expand about the vacuum defined by < X(x) >oc m| ^ 0, where m z is a spinon gap 
(magnon mass), in appropriate units. Then, we can tune the parameter of our system so as 
to define a fully relativistic field theory about the nodes of a c?-wave gap PJ, such that, when 
(3 7^ 0, the following is satisfied: 

< \(x) > (3 - 6/3 + 3/? 2 < FV 3 ^ Q >= (10) 

Note that the non-zero dynamical condensate of the (non- relativistic) operator < \1/ Q cr 3 \l/ Q >, 
obtained above, is compatible with a dynamical opening of a fermion mass gap in the resulting 
relativistic field theory. 

In this way, the fermion terms in the constraint (]5|) decouple, and the effective theory of 
the excitations at the nodes of the ci-wave superconducting gap can be described, up to terms 
that are renormalization-group irrelevant operators in the infrared, by the effective theory (||) 
with Thirring four-fermion interactions, and a standard CP 1 constraint: 

En 2 = i (ii) 

a=l 
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The latter implies that the magnon fields, in their massive (spin gap) phase, can be integrated 
out in a standard fashion in the path integral [pi, to yield an alternative low-energy theory, 



that of a dynamical SU(2) x Us(l) gauge group, with Maxwell kinetic terms for the gauge 
fields, which are the dominant terms in a derivative expansion. 

Superconductivity in this model occurs |I| as a result of dynamical generation of a parity- 
conserving fermion mass in the strong-coupling regime of the Us (I) gauge field fl8| , [5], 0? 



upon coupling the system to external electromagnetic potentials. This dynamical generation 
phenomenon occurs in the infrared region of the effective theory obtained after z-magnon 
integration. In such a theory, upon the opening of a mass gap in the fermion (holon) spectrum, 
the Feynman matrix element: S a =< B®\J u \0 >,a = 1,2,3, with = \&7 /t \I/ the fermion- 
number current, is non-trivial. Due to the colour-group structure, only the massless B 3 gauge 
boson of the SU(2) group, corresponding to the 03 generator in two-component notation, 
contributes to the matrix element. The non-trivial result for S 3 arises from an anomalous 
one-loop graph, depicted in figure 1, and it is given by |25|, |5|: 



S 3 =< Bl\J v \0 >= (sgnm/)^,-^ (12) 



where mj is the parity-conserving fermion mass, generated dynamically by the Us(l) group. 
As with the other Adler-Bell-Jackiw anomalous graphs in gauge theories, the one-loop result 



12) is exact and receives no contributions from higher loops p9 




Figure 1: Anomalous one-loop Feynman matrix element, leading to a Kosterlitz-Thouless-like 
breaking of the electromagnetic Z7b(1) symmetry, and thus superconductivity, once a fermion 
mass gap opens up. The wavy line represents the 577(2) gauge boson B 3 , which remains mass- 
less, while the blob denotes an insertion of the fermion-number current J M = ^7^. Continuous 
lines represent fermions. 



This unconventional symmetry breaking fll2"|), does not have a local order parameter [p5 , [5| 



since the latter is inflicted by strong phase fluctuations, thereby resembling the Kosterlitz- 
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Thouless mode of symmetry breaking^. The massless gauge boson B^ of the unbroken 03 — £7(1) 
subgroup of SU (2) is responsible for the appearance of a massless pole in the electric current- 
current correlator ||, which is the characteristic feature of any superconducting theory. As 
discussed in ref. 0, all the standard properties of a superconductor, such as the Meissner 
effect, infinite conductivity, flux quantization, London action etc. are recovered in such a 
case. The field B^, or rather its dual defined by <9 M = e^ vp d v B z p1 can be identified with the 
Goldstone boson of the broken U em (l) (electromagnetic) symmetry ||. We shall come back to 
the exactness of this result, upon including non-perturbative effects (instantons), in the context 
of our supersymmetric model, later on. 



3 N=l Supersymmetric Gauge Theory of Spin-Charge 
Separation 

We are now ready to discuss the possibility of the emergence of a N = 1 space-time super- 
symmetry in the ansatz ([I]). The main idea behind such a supersymmetrization is to view 
the magnons z as supersymmetric partners of the holons For simplicity, in this note we 
shall turn off the SU(2) interactions in (^), keeping only Us{l), which is mainly responsible for 
the chiral symmetry breaking (mass generation) phenomenon. The incorporation of the gauge 
group SU{2) x UsiX) (@) is straightforward. In this section we shall demonstrate the possibility 
of a N = 1-supersymmetric extension of the action (|J), and of the constraint ((5|), in the absence 
of (non-supersymmetric) external electromagnetic potentials. 

The basic "matter" multiplet of N=l supersymmetry in three space-time dimensions, can 
be written in terms of a scalar superfield as [2E 



® = <j) + 6x + \06F (13) 

which contains a real scalar field, 0, a Majorana spinor \a and a real auxiliary field F. We 
consider complex superfields 

Z= -= ($!+i$ 2 ) =z + 0ip + ^eeF (14) 
v 2 2 

which contain a complex scalar, z = ^75 (0i + ifo), a Dirac spinor, ^ = ^C0i + ^2), and a 
complex auxiliary field, F = ^(-^1 + iF 2 ). The supersymmetry transformations read, 

$sz = & 

S S F = -itffl, (15) 



4 This may be important from a condensed-matter viewpoint, since the absense of a local order parameter 
implies that the opening of a fermion mass gap at the nodes of the original d-wave superconducting gap of the 
cuprate does not affect the d-wave nature of the state. 
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and the supersymmetric invariant lagrangian is given by the highest component (68) of the 
superfield 

DZ*DZ (16) 

where 

Da = W a ~ ima (17) 

is the supersymmetry covariant derivative. 

The gauge field is incorporated in a real spinor superfield which, in the Wess-Zumino gauge, 
takes the form 

v a = iW) a + \ee Va (is) 

where r\ is the supersymmetric partner of the gauge field (gaugino). 

The supersymmetric gauge invariant lagrangian for the matter fields which in terms of super- 
fields is the highest component of the superfield 

VZ*VZ (19) 

with 

V a = D a - iV a (20) 
In terms of component fields the lagrangian reads: 

L = gl\D^z a D^z a + i^pm + F a F a + 2i(f]^ a z a - ^ a rjz a )} (21) 

where denotes the gauge covariant derivative with respect to the Us(l) field, and for con- 
venience we have rescaled the fermion fields \I/ and the auxiliary field F by gi, as compared to 
the non-supersymmetric case, in order to facilitate our superfield formalism. Notice that 
contains a supersymmetric partner of the statistical gauge field Us{l)- From the point of view 
of the ansatz ([!]), this is the defining property of the N — 1 supersymmetric point of the model, 
in the sense that the gauge interaction Us(l) 'doubles' its degrees of freedom. From the point 
of view of the statistical model of ref. I0J, this doubling will only be reflected in the form of 
the effective action, after integrating out the Us(l) field. As explained in ref. 0, this field is 
responsible for yielding fractional statistics to the holons and spinons in three dimensions, and 
as such should be integrated out in the effective action of the physically observable degrees of 
freedom. 

It is important to notice that the constraint @ admits a iV = 1 supersymmetric formulation, 
in terms of the superfields Z a ([14]): 



E = 1 ( 22 ) 

a=l 

Upon integrating out the (non-propagating) a M and gaugino rj fields in a path integral for the 
lagrangian (PT|), and using the constraint (^), it is immediate to obtain the following effective 
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action of holons and spinons in the spin-charge separation ansatz ([!]) at the supersymmetric 
point: 

C s = gXW^d^ + iirp!i a + F a F a + l - (z a d,z a - z a d,z a - itt^**) 2 ] + 

9> \9i 2 E{ + E Zc^Vp + k* a zJF , K + E z a z^ a % + h.c.}} )(23) 



The auxiliary fields F a can be solved by means of the constraint (|22|): 

F a F a = i(^M/ Q ) 2 (24) 



The terms inside the logarithm in ( PB"| ) contain no bare mass terms, but only interaction terms 
among z and $ fields. This can be readily seen by the following formal expansion: 



k=l 



x _ V 2k-l 



\nx = 2 V ~ ~ (25) 



which truncates due to the Grassman structures in x. 

Several important comments are now in order. First, notice that the supersymmetric exten- 
sion of the effective lagrangian for spinons and holons contains both Gross-Neveu and Thirring 
four-fermion interactions. This can be seen by using the Fierz rearrangement formula in three 
space-time dimensions: 

labl^cd = 2<5 a A c - 5 ab 5 cd (26) 
upon which the four fermion Thirring interactions become: 

2 

l^a^a} 2 = -3 E t^A] 2 " 2[^1^2 + 2^2^] (27) 
a=l 

showing that the Gross-Neveu terms in the Thirring interactions cannot cancel the ones ap- 
pearing in ([23|) due to the supersymmetric extension. 

As a result of supersymmetry, the couplings of the four-fermion terms are all related, and 
are of order g\. In the context of the statistical model, such a restriction will imply special 
relations among the microscopic parameters, such as hopping elements, Heisenberg exchange 
energies, doping concentration etc. For instance, in the special case of ref. [0, where the next- 
to- nearest-neighbour (NNN) hopping element t' is assumed dominant, with t ~ 0, one can show 
that four-fermion Gross-Neveu type terms come with generic coefficients of order (t') 2 / (J' S 2 ) , 
with J' (<< J) the NNN Heisenberg exchange energy. In such a situation, supersymmetry 
enforces the relation t' ~ V JJ'5 3 ^ 2 , which may be interpreted as implying that supersymmetric 
points in our formalism may be obtained by tuning the doping concentration. Such restrictions 
may be compatible with the t' = J'/2 supersymmetric point of ref. |13[). In the case above, 
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such an extra restriction implies underdoped situation <5 3 / 2 ~ JJ'/iJ « 1. In more realistic 
models, like the one discussed in ref. ||, involving nearest-neighbour hopping, there will be 
more constraints, involving the hopping element t, etc. A complete analysis along these lines 
falls beyond our present scope. 

Another important issue concerns the physical interpretation of the Majorana fermion rj, 
which, as one can see from ([2iP , fl23D , leads to effective electric-charge violating interactions on 
the spatial planes. From our two-spatial dimension point of view, such violations may admit the 
interpetation of describing interlayer hopping of spin and charge degrees of freedom (hopping 
of real four-space-time-dimensional electrons). In this interpretation, the gaugino rj terms in 
(f21~D, constitute a Majorana-spinor representation of the absence of spin and charge at a site of 
the planar lattice system: 




The reader is advised to draw a comparison with the Grassmann XiX > representation of a 
Wilson line ('missing spin' S ) in the treatment of static holes in refs. [ f27| [J: 

J dx ] dxe~ lS $ ^Eit-^xJxiocKM) (29) 

where a is the temporal component of the gauge potential of the CP 1 er-model, describing 
spin excitations in the antiferromagnet. From this point of view, the existence of N — 1 super- 
symmetry in the doped antiferromagnets necessitates interplanar couplings, through hoping of 
spin and charge degrees of freedom (electrons) across the planes. 



4 Dynamical Mass Generation and N=l Supersymme- 
try 

Next, we proceed to discuss the dynamical scenario for fermion mass generation. First, we 
note that dynamical-mass generation (pairing) in non-supersymmetric models, with combined 
Gross-Neveu and Thirring four-fermion interactions, is possible in three space-time dimen- 
sions. By using a four-component fermion formalism one obtains consistent solutions to the 
Schwinger-Dyson (SD) equations, with non-zero mass, which conserve parity and time-reversal 
invariance p8], ^9| Q. In ref. |29| it was shown that in models with mixed Thirring and Gross- 
Neveu interactions, it is essentially the Gross-Neveu coupling gcN which determines the critical 
behaviour (critical flavour number) of the theory, in a large N expansion. For gcN > 9gni 

5 Note that theories with four-fermion interactions are not in general vector-like, and hence the theorems 
of ref. p0||, for absence of spontaneous violation of parity and time-reversal symmetry due to energetics, 



cannot apply. However, in our superconducting model, integrating out the magnon fields one obtains 24, 23|] a 
dynamical gauge theory in the infrared. It is in this sense that we are interested only in parity-conserving mass 
gaps, which from the point of view of the (low-energy) effective gauge theory, are the energetically preferable 
configurations |3C|. 
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where Qq N is the critical coupling of the (2+l)-dimensional Gross-Neveu model [^T|, the sys- 



tem is dominated by the Gross Neveu interaction, while for gcN < 9gn-> the system becomes 
Thirring like []. 

We now argue that qualitatively the mass -generation phenomenon cannot be affected by the 
presence of supersymmetric partners of the fermion fields. Indeed, the only extra terms in the 
lagrangian (|23| ) that could affect the dynamical mass generation are the terms mixing bosons 
and fermions, zd^z^fj^ . However, at the level of the effective action obtained from (|23"D by 
path-integrating out the z fields, the leading order contributions in a (low-energy) derivative 
expansion, are of order: / d 3 x^/ , ~f fl ^/[(d 2 g fiu — d^d u ) / m 2 z ]^^j v ^ . Such interactions constitute 
irrelevant operators in a renormalization group sense, even at large fermion flavour numbers N, 
and hence do not affect the fixed-point structure of the theory, responsible for mass generation, 
which is thus determined by the four-fermi terms []. 

Within the context of dynamical mass generation, it is important to remark that in super- 
symmetric models dynamical mass generation can occur in a way compatible with unbroken 
supersymmetry only if the effective potential vanishes. This is a result of the equality of the 
fermion and boson masses, mz = rrif = m. In non-supersymmetric theories it is the mini- 
mization of the effective potential that selects the non-trivial solution of the Schwinger-Dyson 
(SD) analysis for the dynamical fermion mass. In contrast, as we shall argue below, in our 
supersymmetric case it is the quantum effective action, and not the effective potential, which is 
responsible for such a selection. The situation is similar to what happens in the two-dimensional 
supersymmetric 0(3) cr-model f3~3fl . In that model, as a result of a constraint similar to (§), 



consistency among the supersymmetry Ward identities, obtained from the quantum effective 
action, selects the non-trivial solution for the dynamical masses, obtained from a SD analy- 
sis [34]. Below we shall not give the details, but we shall present the main arguments, which 
will be sufficient for our purposes in this letter. For simplicity we consider one "complex" 
superfield Z and work with its real components ([111) . The masses of the scalars fa, i = 1, 2 and 
the Majorana spinors x% are related by the Supersymmetry Ward identity: 

< T{ Xl (x),x 3 (0)} >o= (ip + m) < T{fa(x),fa(0) > Q (30) 

where < . . . > Q denote correlators in the non-interacting theory. 
On the other hand, it is known that the fields: 

(-Fi = ~(XjXj), ipXi, d 2 <f>i) (31) 
constitute real superfields T i; the kinetic multiplets of $j. Therefore, the vacuum expectation 



6 In this latter case we should point out that the non-trivial ultraviolet fixed point, found in the numerical 
studies of p2| , might be related - under some sort of ultraviolet-infrared duality - to the non- trivial infrared 
fixed point of the three-dimensional QED, argued in p6[ . 

7 We note that, in a large-flavour-number, N, treatment, these four-fermi operators become renormalizable, 
thereby leading to non-trivial ultraviolet fixed-point structures jn], |8] 
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values of the components of the superfield: 



faTt = (-faFi, XiFi - i<f>ipXi, ' <Pid 2 & + ixM + Wi) (32) 

will be related by the supersymmetry Ward identities. 

Using the equations of motion and the constraint fl22|) this superfield can be written as: 



4>iTi = (u, A, a- z/ 2 ) (33) 

where 

1_ 

A = -i&ftxi 

a = d^d^i (34) 



Then, the corresponding supersymmetry Ward identities become |S3 

S x (p)-(t-2m)D u (p 2 ) = 

j)S x {p) - D a (p 2 ) + 2m{j) - 2m)D u {p 2 ) = (35) 

where D u , D a are the two-point Green's functions of v and a fields respectively, and S\ is the 
corresponding spinorial Green's function of A. Note that the equations of motion, obeyed by 
the Green's functions, have been used in deriving the identities above. 

In the context of the pure Gross-Neveu model in three space-time dimensions, one can 
compute the effective propagators by extending the two-dimensional analysis of ref. [33], in a 



straightforward manner. For instructive purposes we shall derive explicitly the D a propagator, 
pertaining to the Lagrange multipliers a(x) implementing the constraint fl22|) . Expanding about 
the vacuum < a(x) >= m 2 , a(x) =< a(x) > +a'(x), and performing the z integration one 
arrives at an effective action 

S eff>a = J d 3 xTi\n[d 2 + m 2 + a'{x)\ (36) 

The quadratic term in a'(x) determines the effective propagator D a of the quantum field a'. 
Passing onto a Fourier space one obtains: 

S ^ ~JwfJ (Lnk + py-m^'^k^^ (37) 

where m is the dynamically-generated mass for (both) scalars and fermions (due to supersym- 
metry). From this, the propagator D a (p) is obtained immediately. Its p = limit is given 
by: 

D-\0)~f<Pk {h ,* m2)a ~l/m (38) 
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In a similar manner one determines the rest of the Green's functions appearing in (|35|) . For 
the Green's function D u , associated with the linearized Gross-Neveu interactions, we note that 



the quantum corrections have been calculated in ref. |3T|, where a non-trivial ultraviolet fixed 
point structure has been revealed in a large-fermion-flavour number, N, framework. 

With these in mind, one obtains the following results for the pertinent Green's functions, 
to leading order in 1/N expansion: 

m N 
A^O)- 1 ~ m (l + 0[-^ln(A/m)] 



r ii 3 k 1 1 1 

s " (0) ~ / wr w^m^n) + ° 1 n ] ~ non " zero const + (39) 



where A is an ultraviolet cut-off mass. In the above formulas factors of the dimensionful coupling 
constant g\ are understood where appropriate. Moreover, for our purposes in this work the 
detailed form of the 0{jr) corrections will not be important. 
From the Ward identities (1351), on the other hand, one has: 



4m 2 A(0) = -D a {0) 

2mD v (0) = -S x (0) (40) 

Then, on account of (p9|), we see that the first of the identities fl4"0D is satisfied identically to 
this order in 1/N, but one cannot exclude the trivial solution m = 0. Such an exclusion comes 
from the second of the identities (|4"0"D, due to the structure of S\. The so-selected non-trivial 
solution for m, must be the one satisfying the SD equations [j3"H] , by consistency. A non-trivial 
verification of this will come by including the subleading 1/N corrections. The reader should 
keep note of the crucial role of the CP 1 constraint fl2~2"|) in the above selection of the non-trivial 



SD mass gap by the quantum effective action of the N = 1-supersymmetric model |33], p4| . 

In the context of our model, involving both Gross-Neveu and Thirring interactions, a similar 
analysis goes through, with complexities coming from the non-linear realization of supersymme- 
try, and the new interactions in (|2"3"|). Such deviations from the pure Gross-Neveu case, however, 
are in favour of the necessity of a non-zero mass gap, in order to fulfill the supersymmetry Ward 
identities. A detailed analysis along the above lines will be presented in a forthcoming publi- 
cation. For the purposes of this note we restrict oursleves only to pointing out some subtleties, 
associated with the anomalous breaking of the fermion number in our model (c.f. figure 1). In- 
deed, after the z integration, and the implementation of the constraint as above, there are extra 
terms coupling fermions and a multiplier fields in the effective action. One of them involves 
the divergence of the fermion current (after appropriate partial integrations in the action): 

S*f 3 - f d 3 xT A^ r f^ a(x) (41) 



eff J (<9 2 + m 
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If the fermion current Xil^Xi was conserved, then the gauge Ward identity would imply decou- 
pling of this term from the physical correlators. However, as we mentioned above, there are 
anomalies in the model, in the massive phase for the fermions, associated with one-loop graphs 
of figure 1 |5|. |TJ. Such anomalous terms should be properly taken into account in a detailed 
analysis of dynamical mass generation in our supersymmetric model, but we do not expect 
them to affect the selection of the non-trivial solution of the SD equations, characterising the 
pure Gross-Neveu case, studied in detail above. 

An important additional comment concerns the kind of the three-dimensional dynamically- 
generated mass. At present, this seems to depend crucially on the relative sign of the mass, 
between the fermion species. In our analysis above, we have used a single superfield Z, whilst 
in our SU(2) model there are two such superfields. The SD analysis can be extended in that 
case straightforwardly, but alone it cannot make a selection among the two possible signs of 
the mass for these two superfields. Since the four fermion theories are not vector like, one does 
not have at first sight a way of energetically selecting the parity-conserving mass configuration. 
However, as we mentioned previously, the fact that the low-energy integration of magnon fields 
makes the model equivalent to a (vector-like) gauge theory with fermionic matter, is suggestive 



of the exclusion of the parity violating mass, on the basis of the theorems of ref. [30 



5 Instanton Effects, Supersymmetry and Superconduc- 
tivity 

A final issue we would like to address concerns the exactness of superconductivity in the presence 
of non perturbative effects. In the context of the SU(2) x Us(l) theory JTJ, superconductivity 
is associated with the masslessness of the B^ gauge boson of the unbroken U(l) subgroup 
of the SU(2) group, in the massive fermion phase |lj, ||. We now remark that, due to the 
compactness of the pertinent gauge group, instanton configurations - which, in 2+1 dimensions, 
are like monopoles - may give the U(l) gauge boson a small mass. In the dilute- instanton-gas 
approximation, in non-supersymmetric theories, this mass is of order ||35|| : 

m B3 ~ e~ So (42) 

where So is the one instanton action. Such a small mass would destroy the exactness of the 
model's superconductivity, as we remarked earlier. 

We shall argue in this section that supersymmetry favours superconductivity, by further 
suppressing the instanton contributions to the B^ gauge boson mass, as compared to the non- 
supersymmetric case. To this end, we first recall that a dynamical gauge theory is obtained 
in our model by integrating out z and \1/ fields P3J. In a non-supersymmetric theory, upon 
coupling to external electromagnetism, such a procedure leads, in the massive fermion phase, 
to the standard London action for superconductivity |J. In our case, this procedure leads to 
a supersymmetric gauge theory U{1) x Us(l). Indeed, by integrating out z fields one obtains 



Maxwell kinetic terms for the gauge fields, in the phase where the magnon fields are massive |24 . 
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In our supersymmetric theory, the Yukawa coupling of the gaugino 77 to \I/ and z fields results, 
after the integration in Majorana kinetic terms for 77, as required by N = 1 supersymmetry. 



This can be readily seen by a one loop computation, in analogy with the bosonic z part p4 |. 
The relevant graphs, in the massive phase mz = rrif = m (due to supersymmetry), result in 
the following integral: 

^ [ d 3 k K±llL L (43) 

J [{k + p) 2 — m 2 ] [k 2 — m 2 ] 2m 

yielding a Majorana kinetic term ^-rjfir) for the gaugino. One can easily verify the manifest 
N — 1 supersymmetry between this te 
obtained by the z and \1/ integration p4], 



2m 

N = 1 supersymmetry between this term and the corresponding Maxwell terms — ^F 2 U , 



We now remark that in three dimensional supersymmetric gauge theories it is known |26 
that supersymmetry cannot be broken, due to the fact that the Witten index (— 1) F , where 
F is the fermion number, is always non zero. Thus, in supersymmetric theories the presence 
of instantons should give a small mass, if at all, in both the gauge boson and the associated 
gaugino. Although at present there is no rigorous proof of this fact, however, the arguments of 



ref. PB1 indicated that the resulting masses will be even more suppressed than the corresponding 



ones in the non supersymmetric case, 

m Ba = m v = e~ 2So (44) 

with Sn the one-instanton action. 



We should point out, however, that there is an alternative scenario p6 |, which could be 
in operation in our superconducting model. It is possible that supesymmetry is broken by 
having the system in a 'false' vacuum, where the gauge boson remains massless, even in the 
presence of non perturbative configurations, while the gaugino acquires a small mass, through 



non perturbative effects. The life time, however, of this false vacuum is very long [pq| , and hence 
superconductivity can occur, in the sense that the system will remain in that false vacuum for 
a very long period of time, longer than any other time scale in the problem. 

Whichever of the two scenaria is realized in the model, from a condesed-matter point of 
view the important conclusion, obtained from the above analysis, is that the coupling of the 
superconducting planes due to interlayer electron hopping, associated with the presence of the 
gaugino field r], helps stabilizing superconductivity, which otherwise would be jeopardized by 
non-perturbative effects. 



6 Discussion 

In this work we have demonstrated the possibility of N=l Supersymmetric gauge theories in 
the context of a spin-charge separation ansatz of the SU(2) x Us(l) gauge model of ]T]. Such 
models may be relevant for the physics of superconducting gaps which open up at the nodes of 
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a d-wave gap of high-T c cuprates. The super symmetry was achieved without introducing un- 
physical degrees of freedom. However it necesitates the coupling of the superconducting planes. 
Its presence, seems to suppress the effects of instantons of the gauge field SU(2), which could 
jeopardise superconductivity in the model. As far as the lattice system is concerned the super- 
symmetry is achieved modulo irrelevant operators in a renormalization-group sense. This may 
imply that our considerations in this work might also be relevant to the construction of more 
general supersymmetric gauge theories on the lattice, in the sense of obtaining supersymmetric 
continuum theories by droping possibly non-supersymmetric, renormalization-group irrelevant 
operators. 

We believe that our work may prove useful towards an exact discussion of phase diagrams of 
three-dimensional effective gauge models of antiferromagnetic superconductors, via the analysis 
of the quantum moduli space of gauge theories, in the spirit of Seiberg and Witten fp~4]| . In this 
respect, we note that N — 1 supersymmetry in three dimensions, which we have considered here 
as the minimal way of supersymmetrization of a doped spin-charge separated antiferromagnet, 
without the introduction of extra degrees of freedom, cannot yield exact results. It is the 
N = 4 supersymmetry in three-dimensional gauge theories which can produce such results. In 
three dimensions, N = 2 theories may also allow for some exact results, in connection with the 
geometry of their quantum moduli space |p"5 |. 

At present, our physical understanding for a condensed-matter spin-charge separated model 
exhibiting N = 2 supersymmetry is not complete. One might speculate that, since N = 2 
three-dimensional supersymmetric theories are obtained |H| by dimensional reduction of four- 
dimensional N — 1 supersymmetric theories, such models might have some relevance to a 
possible extension of the ideas in our work beyond the planar structures. We shall present a 
more detailed study of such models in a future work f36[. We should remark however that, 



as far as the spinon and holon degrees of freedom are concerned, the extension to N = 2 
supersymmetry is immediate, with no extra doubling of degrees of freedom. The novel feature, 
compared to the N=l case, is the presence of a Dirac-like gaugino. Due to its Dirac nature, 
the gaugino may now carry non-trivial charge under the external electromagnetism, and thus 
the effective action conserves the electric charge, in contrast to the present situation with a 
Majorana gaugino. In view of the aforementioned embedding of 3-dimensional N=2 theories 
in 4-dimensional N=l supersymmetric theories, this possibility of conservation of the electric 
charge may be related to the exact conservation of electric charge in four dimensional space 
times. From the point of view of dynamical mass generation, we should remark that, at first 
sight, the N = 2 3-dimensional models appear not to generate a dynamical mass. This is due to 
the fact that such theories are obtained from N = 1 4-dimensional models, where claims have 
been made |37| that non-renormalization theorems in the supersymmetric SD equations yield 



only the trivial solution for the mass. However, such claims have been questioned recently p8 
From our point of view we consider the issue as still open. 

Another comment we would like to make concerns the fate of supersymmetry at finite 
temperatures. We expect the supersymmetry to be broken at finite temperatures, which results 
in different masses for spinon and holons, a situation probably met in realistic cases. However, 
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even in such a case of broken supersymmetry, the existence of a supersymmetric vacuum at 
zero temperatures is useful in providing some exact information about the phase diagram along 
the lines mentioned above. 

Before closing we should also stress that our results apply even to one- dimensional chains of 



holons, which may characterize certain underdoped cuprates in the so-called stripe phase [39 
Such systems appear to be described by a spin-charge separated phase, where the holon de- 
grees of freedom lie on one-space dimensional stripes (chains), spatially separated by regions 



of zero doping. As discussed in ref. |27|], spin charge separation in one (spatial) dimensional 
antiferromagnetic models leads to gauge theories of Dirac fermions coupled to a CP 1 a- model. 
The continuum action is similar in form, but in two space-time dimensions, with the action 
([|). In such a case, the resulting N — 1 supersymmetric extension will again involve combined 
Gross-Neveu and Thirring interactions. Such (l+l)-dimensional models have been studied pre- 



viously in the literature [S3]. As far as supersymmetry and dynamical mass generation are 
concerned, such models share the same qualitative features as their (2 + l)-dimensional coun- 
terparts, discussed here. The gauginos in such one-dimensional theories could then describe 
(effective) electron hopping across the chains. At present we are agnostic as to whether such 
supersymmetric spin-charge separated models play any crucial role on the physics of the stripe 
phase of the underdoped cuprates. 
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